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ABSTRACT 

W28 is one of the archetype supernova remnants (SNRs) interacting with molecular 
clouds. H.E.S.S. observation found four TeV sources which are coincident with the 
molecular clouds (MCs) around W28, but Fermi LAT detected no prominent GeV 
counterparts for two of them. An accumulative diffusion model is established in this 
Letter and the energetic protons colliding the nearby MCs are considered to be an 
accumulation of the diffusive protons escaping from the shock front throughout the 
history of the SNR expansion. We have fitted the 7-ray spectra of the four sources 
and naturally explained the GeV spectral break of the northeastern source (source N) 
and the nonsignificant GeV emission of the southern sources A and C. The distances 
of sources A and C from the SNR centre are found to be much larger than those of 
sources N and B, which may be the basic reason for the faint GeV 7-rays of the two 
former sources. 

Key words: radiation mechanisms: non-thermal 
pernova remnants. 
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1 INTRODUCTION 

Cosmic rays (CRs) below the "knee" in the Galaxy are com- 
monly believed to be accelerated at the shock of supernova 
remnants (SNRs). However, whether the 7-rays from SNRs 
are of hadronic or leptonic origin is still in hot debate. Al- 
though evidences for 7-ray emission arising from molecular 
clouds (MCs) have been found in two SNRs, IC 443 (Ac- 
ciari et al. 2009) and W28 (Aharonian et al. 2008), it is 
not easy yet to confirm the emission of the accelerated pro- 
tons, since the competing lepton processes can also account 
for the high energy radiation. Fortunately, emission from 
MCs illuminated by CRs from nearby SNRs (Aharonian & 
Atoyan 1996; Gabici et al. 2009) can give us opportunities 
not only to distinguish the contribution of the hadrons from 
that of leptons but also to investigate the escape process 
of these energetic protons (Fujita et al. 2009). SNR W28, 
with recent observation in GeV and TeV 7-rays, is such an 
excellent case. 

W28 (G6.4— 0.1) is one of the prototype thermal 
composite (or mixed-morphology) SNRs, characterized by 
centre-filled thermal X-ray and shell-like radio emission, 
with large (about 50' x 48') apparent dimensions (Green 
2009). It is considered to be an evolved SNR in the radia- 
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tive stage, with an age estimate spanning from 35 to 150 kyr 
(Kaspi et al. 1993). SNR W28 has been confirmed to be in- 
teracting with MCs with robust evidences including the 1720 
MHz OH masers, morphological agreement with molecular 
features, molecular line broadening, etc. (see Jiang et al. 
2010; and references therein). The observations of molecu- 
lar lines place it at a distance of ~ 2kpc (Velazquez et al. 
2002). 



In 7-rays, H.E.S.S. has revealed four TeV sources in 
the W28 field positionally coincident with MCs (Aharonian 
et al. 2008): HESS J1801-233 (denoted as source N in this 
Letter), located along the northeastern boundary of W28, 
and HESS J1800-240A, B, and C (denoted as sources A, B, 
and C, respectively, in below), located ~ 30' south of W28. 
Recently, Fermi LAT observation (Abdo et al. 2010) found 
GeV sources 1FGL J1801.3-2322c and 1FGL J1800.5-2359c, 
which coincide with sources N and B, respectively. More 
interestingly, it did not detect significant GeV counterparts 
for sources A and C. The LAT upper limits combined with 
the H.E.S.S. data imply a spectral break between lOGeV 
and lOOGeV. We will show in this Letter that these 7-ray 
spatial and spectral properties can be accounted for with an 
accumulative diffusion process of CRs which are accelerated 
at the dynamically evolving blast shock front. 
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Figure X. Sketch for the positional relation between the source 
point (S) and the field point (C) (see text in 



2 ACCUMULATIVE DIFFUSION MODEL 

CR diffusion models have been developed by Aharonian & 
Atoyan (1996) and Gabici et al. (2009) to evaluate the CR 
spectrum in the vicinity of SNRs. In their models, both 
cases of impulsive and continuous injection are considered 
on the assumption that the distance between the SNR and 
the cloud is significantly larger than the size of the SNR, so 
that the injection region can be approximated as a point- 
like source. If, however, the cloud is at or close to the SNR 
blast shock, such as in the case of W28, the diffusion dis- 
tance for protons injected at various time, ti, and various 
shock radius, R s , is not a constant and the assumption of 
point-like source injection should be modified. Therefore, we 
here establish an accumulative diffusion model considering 
that the energetic protons colliding the given MC are a col- 
lection of the diffusive protons escaping from different _R a as 
the SNR expands. Small distance between the SNR and MC 
is allowed. 

The distribution function at an arbitrary field point C 
(at radius R c ) of the energetic protons that escape from unit 
area at an arbitrary source point S on the spherical shock 
front surface (see Figure[T]) is given by (Aharonian & Atoyan 
1996) 
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where Rfe,0,<j>) = \/Ri(ti) + R% - 2# s i? c sin 9 cos <f> is the 
distance between points S and C (Figure [TJ , tdif = i ag c — ti 
is the diffusion time after the escape, Rdn = Rdif(E p , tdu) = 
2-v/ D(E p )tdit is the diffusion radius, and dNo/dti = QoE~ F 
is the total escaping rate at time ti. Here the diffusion 
coefficient is assumed to be in the form of D(E P ) — 
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where x is the correction fac- 



tor of slow diffusion around the SNR (Fujita et al. 2009) 
and 8 « 0.3-0.7 (Berezinskii 1990). The number of pro- 
tons with energy E p escaping from the SNR shock front 
surface at ti and arriving at the spherical surface of radius 
Rc at t agc is A-kR 2 s J^Q(E p ,R(ti, 9, cf>), t d n) R 2 sin 



Because of spherical symmetry, the distribution func- 
tion of protons (escaping from Rs(ti)) for any point 
(e.g., point C) on the spherical surface of radius R c is 
R 2 J^£f(E p ,R(ti,e,(j}), tdu) sin 9d9d(j>. Therefore, the dis- 
tribution function for point C accumulating all the historical 
contributions can be obtained as 
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For the dynamical evolution of radius of the SNR ex- 
panding in the interstellar (intercloud) medium of den- 
sity po — l.4m H »io, we use the Sedov- Taylor law R s — 
(2. 026 E SNR t 2 1 po) 1 ' 5 for the adiabatic phase and _R S = 
(14:7eE SNR R 2 t 2 /inpo) 11 " 7 for the radiative phase, where -Rt 
is the transition radius from the Sedov phase to the radiative 
phase, -Esnr is the supernova explosion energy, and e is a fac- 
tor equal to 0.24 (Blinnikov et al. 1982; see also Lozinskaya 
1992). 

The total energy of the escaping protons is W p = 
vEsnr = tage J QqE~ v E p dE p , where r\ is the fraction of 
the explosion energy converted into protons. 

In addition to the escaping protons, we also consider the 
contribution of diffuse Galactic protons in the Solar neigh- 
borhood with flux density (see, e.g., Dermer 1986): 



Jcr(£ p ) = 2.2 
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In the calculation of the 7-rays from the nearby MC (of 
mass M c ) due to p-p interaction so as to match the observed 
GeV and TeV fluxes, we use the analytic photon emissivity 
dN 1 /dE 1 developed by Kelner et al. (2006). At high ener- 
gies, 

dNj 
dE~, 
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where rib is the average density of target baryons of the 
cloud, Cinei is the cross section of inelastic p-p interactions 
and function F 7 (x, E p ) is defined by Eq. (58) in Kelner et al. 
(2006); while at low energies, with the modified 5-function 
approximation as proposed in Aharonian & Atoyan (2000), 
the emissivity of 7-rays is 
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where i? m i n 
mesons is 



E-y + m 2 /4Ej, and the production rate of 7r° 
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(see Kelner et al. 2006 for the explanations of other param- 
eters). 

The contribution of 7-rays from the secondary leptons 
that are created by p-p interaction is negligible, as compared 
with the dominant contribution of the protons themselves 
(Gabici et al. 2007; 2009). Because the timescale for the sec- 
ondaries to escape from the cloud is shorter than the energy 
loss time for particle energies between ~ 100 MeV and a 
few 100 TeV, a large amount of the secondaries produced 
within clouds of a wide range of parameters escape without 
being affected by significant losses (Gabici et al. 2009). In 
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Table 1. Parameters of the model considered for different sources. 



Variable 


i? c /pc 


8 


M C /10 4 M Q 


nt/cm 3 a 


source N 


12 


0.45 


4 


1.6 x 10 3 


source B 


22 


0.35 


0.2 


1.1 x 10 2 


source A m i n 


100 


0.6 


3.2 


1.8 x 10 3 


source A 


140 


0.65 


8 


4.4 x 10 3 


source A max 


180 


0.7 


13 


7.2 x 10 3 


source C m i n 


55 


0.55 


0.64 


1.0 x 10 3 


source C 


100 


0.6 


2.4 


3.9 x 10 3 


source C max 


250 


0.7 


22 


3.5 x 10 4 



a The adoption of the average angular diameter of each source 
for deriving the baryon density is described in the text in [[3] 



the molecular gas of density as high as 10 4 cm -3 , the mag- 
netic field could be assumed to be of order 10 2 /iG (Crutcher 
1999). Extremely high energy (> a few 100 TeV) secondaries 
radiate all their energy by synchrotron peaked at tens of keV 
before escaping from the cloud; but the X-rays are not in the 
scope of this letter. 



3 APPLICATION TO SNR W28 

Now we use the accumulative diffusion model established 
above to reproduce the GeV-TeV spectra for the four 7-ray 
sources around SNR W28, showing that the different 7-ray 
spectra of the four sources mainly result from the different 
spectra of the runaway protons, as a consequence of the 
different distances of the sources (molecular clouds) from 
the SNR shock from which the accelerated protons escape. 

The current radius of W28 is approximated as 
Rsnr ~ 12 pc by using 2kpc as the distance to it. 
The current blastwave velocity usnr was inferred to 
be in the 60-100 km s -1 range based on the optical 
observations (see Rho & Borkowski 2002), and hence 
we adopt usnr = 80kms _1 . Also following Rho & 
Borkowski (2002), we adopt the intercloud medium den- 
sity n ~ 10 cm -3 . For the evolution in the radiative 
phase, its age is estimated as i agc ~ 2R SNR / (7v S nr) = 
4.2 x 10 4 ( J R SN R/12pc)(7; SN R/80kms~ 1 )- :L yr- The explo- 
sion energy is given by -Esnr = 6.6 x 10 5O (n o /10 cm 3 ) 
(«SNR/80kms~ 1 ) 1,2 (_R S NR/12pc) 3 erg. Moreover, we assume 
a power-law energy spectrum of the escaping protons of in- 
dex T = 2.2 (e.g., Giuliani et al. 2010), and adopt x = 0.1 
(Fujita et al. 2009) and rj = 0.1 (Blandford & Eichler 1987). 
Therefore, there are three adjustable parameters: R c , S, and 
M c . These model parameters for the 7-ray sources around 
SNR W28 are shown in Table 1. All the spectra of the four 
sources are well reproduced (Figure[2]), as described in detail 
below. 

Source N is a well known shock-MC interaction re- 
gion on the northeastern boundary of SNR W28. We adopt 
i? c ~ -Rs(tagc) ~ 12pc for this source. The observed GeV 
spectrum of source N has a spectral break at 1.0 ± 0.2 GeV 
with photon index 2.1 below the break and 2.7 above the 
break (Abdo et al. 2010). As seen in the top left panel of 
Figure [2] this spectral break together with the TeV spec- 
trum can be naturally explained by the accumulative dif- 
fusion model without invoking either a breaking power-law 



proton spectrum (see Abdo et al. 2010) or any artificial low 
energy cut-off. The model spectrum peaks at ~ 0.35 GeV. 
We note that the TeV emission with a slope 2.7 implies the 
same power-law index for the protons, and this is unlikely to 
come directly from the shock accelerated protons (because 
the typical accelerated particle spectral slope for accelera- 
tion process is 2-2.3) but can be obtained from a diffusion 
process. 

Source B is located right on the south of SNR W28 
with a projected distance of ~ 10 pc from its southern 
circular boundary (Aharonian et al. 2008), and hence we 
adopt R c ~22 pc. In source B and also in source N, which 
are closely near the SNR, the contribution from the diffuse 
Galactic CRs is negligible compared with that of the diffused 
protons from the SNR. In our model, with the increase of 
radius R c and hence the increase of the mean diffusion dis- 
tance, the spectral peak shifts to higher energies. The peak 
~ 0.62 GeV for source B (see the top right panel of Figure(2]| 
is higher than the peak in source N, which implies the diffu- 
sion effect for different diffusion distances. Source B is also 
suggested to be associated with HII region W28A2 which 
may also possibly explain the 7-ray emission, but requires 
extremely high density of gas (Abdo et al. 2010). 

Sources A and C were not found to have significant 
counterparts in the Fermi LAT observation. The LAT up- 
per limits together with the H.E.S.S. data are indicative of a 
spectral break in the 10-100 GeV range. This break is repro- 
duced in Figure [2] (bottom left panel), in which the concave 
7-ray spectra appears as "M" shape with two peaks. The 
underlying proton spectra consist of two components: the 
galactic proton background, which is responsible for the left 
peak, and the protons coming from SNR W28, which are 
for the right one. With the increase of R c , the flux density 
of CRs accelerated by the SNR is diluted and, compared 
with the SNR protons, the relative significance of the dif- 
fuse Galactic protons rises. There are uncertainties in the 
GeV fluxes for both sources A and C. The lower limits of R c 
(dotted lines in the bottom left panel of Figure [2]) are used 
for the 7-rays produced by protons escaped from SNR to 
match the upper limits of the LAT data, and the upper lim- 
its of -Rc (dashed line in the bottom left panel of Figure [2| 
are for the diffuse Galactic protons to match the LAT up- 
per limits. Meanwhile, with R c increased, the (right) peak 
of the emission of the escaping SNR protons now shifts to 
higher energy than that for source B, thus forming the spec- 
tral break at 10-100GeV. Although there is less confinement 
in the GeV realm for source C, we can determine the up- 
per and lower limits of R c of the source (see the dotted and 
dashed lines in the bottom right panel of Figure 01 , in a 
similar way used for source A. The fitted values of R c for 
the both sources are much larger than those of sources N 
and B (see Table 1), which maybe the basic reason why the 
GeV emission for A and C looks faint. 

The estimated mass of northeastern cloud (correspond- 
ing to the 7-ray source N) from NANTEN data is ~ 5 x 
10 4 Mq and the total mass of the southern clouds (includ- 
ing the clumps corresponding to sources A, B, and C) is 
~ 1O 5 M (Aharonian et al. 2008). The values of M c fitted 
for the four sources in our model (Table 1) seem to be in 
agreement with the NANTEN observation. Using the aver- 
age angular diameters of sources N, A, B, and C, 18', 16', 
16', and 12', respectively, as measured from the outermost 
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Figure 2. 7-ray spectral data of the four sources and the model spectra. The observed GeV and TeV data are adapted from Abdo et 
al. (2010) and Aharonian et al. (2008), respectively. The model parameters are given in Table 1. The dashed and dotted lines represent 
the model spectra for the upper and lower limits of R c , respectively, for both sources A and C (see text in . 



contours of the TeV emission in Figure 1 in Aharonian et al. 4 CONCLUSION 
(2008) , we give the baryon density for each molecular cloud 
in Table 1 also. 



The reproduction of the four 7-ray spectra produces 
various values of parameter 8 in the range of 0.35-0.7. The 
diffusion process of CRs results from particle scattering on 
random MHD waves and discontinuities. Theoretically, the 
power-law index of diffusion, S, then depends strongly on 
the spectral energy density of interstellar turbulence which 
would originate from fluid instabilities and magnetic field 
fluctuations. However, the mechanisms of the turbulent phe- 
nomena are still in hot debate. In general, three types of 
isotropic diffusion scalings (Kolmogorov, 5 = 1/3; Kraich- 
nan, 5 = 1/2; and Bohm, 5 = 1) may exist in the ISM (see 
Strong et al. 2007). The S values obtained in our model vary 
just between these values and may arise from one or a com- 
bination of the three types of diffusion, in view of the fact 
that the escaping protons diffuse to the four MCs though a 
nonuniform and anisotropic interstellar environment around 
SNR W28. 



An accumulative diffusion model has been established and 
well explains the 7-ray spectra of the four sources around 
SNR W28 recently obtained by H.E.S.S. and Fermi. The 
energetic protons colliding the nearby MCs are considered 
to be an accumulation of the diffusive protons escaping from 
the shock front throughout the history of the SNR expan- 
sion. For the various distances of the MCs from the SNR 
centre, the resulted proton spectra can have prominently 
different shapes. We have thus fitted the 7-ray spectra of 
the four sources and naturally explained the GeV spectral 
break of the northeastern source N and the nonsignificant 
GeV emissions of the southern sources A and C. The dis- 
tances of sources A and C from the SNR centre are found 
to be much larger than those of sources N and B. This also 
implies prospects for inferring the spatial position of the 
proton-illuminated MCs around SNRs using the GeV and 
TeV spectra. 7-rays from molecular clouds illuminated by 
accumulated diffusive protons from SNR W28 strongly sup- 
port the scenario that hadrons accelerated by SNR shocks 
contribute to the Galactic CRs. 
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